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We describe how the Schrodinger functional and twisted mass QCD can be used to compute the kaon B- 
parameter Bk on the lattice with Wilson fermions. This new approach is expected to reduce the systematic 
uncertainties on Bk through a better control of its (nonperturbative) renormalisation. Preliminary results for 
the bare matrix element in a physical volume of 1.5"^ x 4.5 fm* are presented. The renormalisation of the matrix 
element in a Schrodinger functional scheme is also discussed. 



1. INTRODUCTION 

Out of the physical quantities related to CP 
violation in the Standard Model, the parame- 
ter e accounting for indirect CP violation in the 
kaon system is the one with the simplest non- 
perturbative contribution. It is given by the kaon 
_B-parameter Bk, defined as: 

- W^^2 ' 

where O'^^^^ is the AS" = 2 operator of the ef- 
fective weak Hamiltonian: 

OA5=2 ^ ^-,^Ld){-S^Ld) , (2) 

and 7^ — ^7/i(l — 75). In an obvious notation, it 
is convenient to split O^^^^ in parity-even and 
parity-odd parts: 

^VA+AV ■ (3) 

Since parity is a QCD symmetry, only the parity- 
even part O0y^\j^ will yield a nonvanishing K^- 
matrix element. 
In spite of the absence of both power-divergent 
subtractions and final state interactions, the lat- 
tice computation of Bk is still a very difficult 

*Based on a talk by C. Pena at the International Sympo- 
sium of Lattice Field Theory, August 19-24 2001, Berlin, 
Germany. 
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problem. Besides the systematic error arising 
from the quenched approximation, the breaking 
of chirality (due to the use of Wilson fermions) 
leads to a mixing under renormalisation of Q^^^^ 
with four other operators of dimension six. It 
is now understood that the corresponding mix- 
ing coefficients, as well as the overall, scale- 
dependent, renormalisation constant, require an 
accurate nonperturbative computation. This has 
been dealt with by using Ward identities and RI- 
MOM renormalisation |l|,|| , but the result still re- 
flects large numerical uncertainties. On the other 
hand, the use of lattice actions with better chi- 
ral properties has yielded results which have not 
manifested clear agreement with those obtained 
with Wilson fermions. We refer the reader to ||] 
for reviews on the subject and references. 

To improve on this situation, we propose a new 
computation of Bk with Wilson fermions. It 
adopts the Schrodinger functional (SF) Q for the 
computation of scale-dependent renormalisation 
constants, and a twisted mass QCD (tniQCD) ac- 
tion iQ, in order to avoid the mixing of O^^^^ 
with other operators under renormalisation. An 
alternative proposal for the elimination of this 
mixing has been presented in 
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2. tmQCD FOR Bk 

Twisted mass QCD was proposed as a means 
of avoiding the problem of exceptional configu- 
rations, allowing lattice QCD computations with 
light quark masses closer to the chiral limit H]. 
For the computation of Bk we will use an action 
with a twisted doublet of light quarks "0 — {u, d) 
and a standard (untwisted) strange quark s: 

S = a'^Y^ [{p{x) {D + mi+ mj5T^) Tpix) 

X 

+ s{x){D + ms)s{x)\ , (4) 

where acts on isospin indices and /i; is 
the "twisted mass". In the classical theory 
at the continuum limit, the standard QCD ac- 
tion is recovered by changing fcrmionic variables 
through an axial rotation of the form ?/; — > 
exp{— ia75T'^/2}V', where tan(a) = ^i/mi. This 
rotation leaves invariant the combination Mi — 
yjmf + ^f, which can therefore be regarded as 
the physical bare light quark mass; on the other 
hand, it induces a mapping between composite 
operators in the two theories. It has been shown 
in Q that this classical equivalence carries over to 
appropriately renormalised correlation functions 
in the quantum theory. The twisting angle is then 
given by tan(Q!) = Mi.-R/'^i.i?; where the masses 
are renormalised. 

By fixing renormalised masses so that a — 7t/2, 
the mapping between operators induced by the 
corresponding axial rotation leads to the relation: 



{K"\0 
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0\ct=7r/2 _ 
/tmQCD ~ 



(5) 



(the above is to be understood as valid for renor- 
malised matrix elements). Parity violation at 
nonzero twisted mass ensures that the tmQCD 
matrix element does not vanish. Computing the 
Ihs of Eq. (H), rather than the rhs, is advanta- 
geous, as the parity-odd operator Oy^+^y is pro- 
tected by CVS symmetry from mixing with other 
operators, and hence it renormalises in a purely 
multiplicative fashion This is clearly true if 
renormalisation is carried out in the chiral limit. 



where tmQCD and QCD (and their symmetries) 
coincide. 

We show the feasibility of a direct computa- 
tion of the matrix element of Oyj^^^y with the 
tmQCD action and SF boundary conditions by 
presenting preliminary results for the bare value 
of Bk, obtained on a 16^ x 48 lattice at /3 = 6.0. 
This corresponds approximately to a physical vol- 
ume of 1.5"^ X 4.5 fm''. In the computation we have 
kept the s and d masses degenerate and close to 
the physical s quark mass, as done in previous 
computations of Bk- We show in Fig. ^the ratio 
of correlations: 



-^{0'afirAfAv{x)Ods) 
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(6) 



where IC = —-^570(1 + i"f^)d is the operator ob- 
tained by axially rotating the QCD Aq current 
with a = 7r/2, and 0,0' are the standard SF 
pseudoscalar sources: 



(7) 
(8) 



where /i, /2 label two distinct flavours. The cen- 
tral plateau of R corresponds to the bare value 
of Bk- As expected, it appears in the time in- 
terval 0.35 < xo/T < 0.65, where the single kaon 
states coming from Euclidean times and T are 
isolated (as shown by the effective mass plots of 
Fig. 1). Fitting the i?-plateau to a constant we 
get BK{a-^ ~ 2 GeV) = 0.96(2). 

3. SF RENORMALISATION 

The renormalisation constant of the matrix el- 
ement at renormalisation scale /i — can 
be evaluated on symmetric lattices of physical 
volume L"*. The SF renormalisation condition, 
which is imposed at the chiral limit m/./f = 
rUs^B. = ^J■l,B. = 0, is then: 
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16=x4B, ;S = 6.0, (k,,kJ = (0. 1342,0. 13518), au =0.0266 
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Figure 1. Plateaux for Bk and effective kaon 
masses, from 250 configurations, for a 16^ x 48 
lattice at /? ~ 6.0. 



where /i is a wave function renormalisation fac- 
tor, defined as: 



(10) 



and F is a correlator involving Oy^+^iy, chosen 
so as not to vanish in the chiral limit. We opt for 
the simple choice: 

Fixo) = L^O[,,0^i=;\y{x)OMu) ■ (11) 

Once ZvA+AV has been defined at given values 
of the bare coupling go and of the renormalisa- 
tion scale n = , a step scaling function that 
describes the (nonperturbative) running with /i 
can be defined at fixed value u = g'^{L) of the 
renormalised coupling as: 

ZvA+Avigoi 2L) 



f^VA+Aviu) = lini 



Z' 
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(12) 



g^(L)=u 



Having computed cryA+AV for a number of values 
of M, it is then possible to run the renormalised 
matrix element to large scales /x = L~^, at which 
the perturbative matching to other renormalisa- 
tion schemes, such as MS, can be performed. 

By computing various Z constants, obtained 
from different sources O, O' (i.e. different renor- 
malisation conditions), wc hope to get an idea of 



the size of 0{a) effects, which arise due to the use 
of the unimproved four-quark operator. 

4. CONCLUSIONS 

We have described how the SF and tmQCD can 
be used to achieve a better control of the renor- 
malisation properties of the kaon i?-parameter 
Bk- a first feasibility study for the bare matrix 
element giving Bk in tmQCD has been presented, 
and the strategy for the SF renormalisation dis- 
cussed. 

We mention that the same techniques can be 
extended to all AS" = 2 operators of interest, thus 
avoiding the spurious mixing with operators of 
the same dimension induced by the use of Wilson 
fermions. The computation of the step scaling 
functions for the complete basis of such operators 
in already under way. 
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